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X-RAY EMISSION FROM THE TYPE IC SUPERNOVA 19941 OBSERVED WITH CHANDRA 



(N 
O 
O 
(N 



(N 



> 

m 

t> 
m 
in ■ 

o : 

(N . 
O ■ 

^ ' 

Or 

6 ■ 

a: 



Stefan Immler 

Astronomy Department, University of Massachusetts, Amherst, MA 01003 

Andrew S. Wilson & Yuichi Terashima 1 

Astronomy Department, University of Maryland, College Park, MD 20742 
Draft version February 1, 2008 

ABSTRACT 

We present two high-resolution Chandra, X-ray observations of supernova (SN) 19941 which show, for 
the first time, that the interaction of the blast wave from a Type Ic SN with its surrounding circumstellar 
material (CSM) can give rise to soft X-ray emission. Given a 0.3-2 keV band X-ray luminosity of 
L x ~ 1 x 10 37 ergs s" 1 between six and seven years after the outburst of SN 19941, and assuming 
the X-ray emission arises from the shock-heated CSM, we derive a pre-SN mass-loss rate of M ~ 1 x 
1CT 5 M Q yr- 1 (u w /10 km s" 1 ). Comb ining the results with earlier ROSAT observations, we construct 
the X-ray lightcurve of SN 19941. A best-fit X-ray rate of decline of L K oc t~ s with index s ~ 1 and 
a CSM density profile of p csm oc r^ 1 - 9 * - 1 are inferred, consistent with what is expected for a constant 
mass-loss rate and constant wind velocity profile for the SN progenitor (p csm oc r~ 2 ). 



Subject headings: supernovae: individual (SN 19941) 
19941) — X-rays: ISM 

1. INTRODUCTION 

The classification scheme for SNe is based on the pres- 
ence or absence of hydrogen lines in their optical spectra 
(Type II and Type I, respectively). While Type Il/Ib SNe 
occur during the core collapse of a massive star, Type la 
SNe are believed to originate from the deflagration or deto- 
nation of an accreting white dwarf star in a binary system. 
The massive (> 10M Q ) progenitor stars for the controver- 
sial Type Ic SN subclass are thought to have lost their 
outer hydrogen and helium layers either through stellar 
winds or by mass transfer to a companion, leaving behind 
a stripped carbon and oxygen (C+O) star. Such a 'naked' 
C+O star will later explode when its iron core collapses 
(e.g. Nomoto et al. 1994). The interaction of the out- 
going SN shock wave with the ambient CSM, deposited 
either by a pre-SN stellar wind or non-conservative mass 
transfer to a companion, produces hot gas with character- 
istic temperatures in the range T ~ 10 7 -10 9 K (Chevalier 
& Fransson 1994). Gas heated to such high temperatures 
produces radiation predominantly in the X-ray range. X- 
ray emission from this interaction is expected for all Type 
Ib/c and II SNe with substantial CSM established by the 
massive progenitors. Over the last 20 years, searches for 
X-ray emitting SNe have been successful for only a rel- 
atively small number of Type II's in the near aftermath 
(days to months) of the explosion (see Immler & Lewin 
2002 for a review article) 2 . 

By contrast, no Type la or Ib/c SN has ever been firmly 
detected in X-rays, apart from a 3.5ct excess in an XMM- 
Newton EPIC-PN image 4" offset from the optical posi- 
tion of the Ic SN 2002ap (Rodriguez-Pascual et al. 2002). 
Hard (2-10 keV band) X-ray emission was recorded with 
BcppoSax from the position of the unusual SN 1998bw, 
which might be associated with a 7-ray burst event (GRB 
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980425; Galama et al. 1998, Pian et al. 1998). Due to 
the large error box of the BeppoSax Wide Field Camera 
observation (between 3' and 8', 99% confidence limit) and 
the rather large probability of ~ 60% that the source is a 
random chance coincidence, the association of the X-ray 
source with SN 1998bw and GRB 980425 is still tentative 
(Galama et al. 1998). 

Evidence for soft (0.1-2.4 keV band) X-ray emission 
has been reported from the Type Ic SN 19941, based 
on ROSAT HRI observations 82 days after the outburst 
(Immler, Pietsch & Aschenbach 1998). However, the 
ROSAT observations, with a spatial resolution of ~ 5" 
(FWHM on-axis), were not conclusive since SN 19941 is lo- 
cated close to the X-ray bright nucleus of the host galaxy 
M51 (distance ~ 18") and is embedded in a high level 
of extended X-ray emission from hot gas and unresolved 
point-like X-ray sources in the bulge of M51. 

2. X-RAY OBSERVATIONS AND ANALYSIS 

We used the Chandra X-ray observatory to search for X- 
ray emission from SN 19941 in X-ray images taken between 
six (June 6, 2000) and seven years (June 23, 2001) after 
the explosion (March 31, 1994; Chandler, Phillips & Ru- 
pen 1994). Our 14.9 ks and 26.8 ks Chandra observations 
were carried out with the S3 chip of the Advanced CCD 
Imaging Spectrometer (ACIS-S) at the focal plane of the 
telescope and the nucleus of the host galaxy M51 placed 
at the aim-point (Terashima & Wilson 2001). The superb 
spatial resolution (0'.'5 FWHM on-axis), together with the 
high sensitivity of the instrument and the long combined 
exposure (41.7 ks), allows us to separate point-like sources 
from the diffuse emission in the bulge of M51 and to carry 
out a sensitive search for X-ray emission from SN 19941. 
The co-added and adaptively smoothed Chandra image of 



the central region of M51 is presented in Fig. 1. 
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An X-ray source is detected at the position of SN 19941 
in the merged Chandra data with a significance of ~ 
6a in the 0.3-2 keV band using the 'wave detect' algo- 
rithm implemented in the CIAO data analysis package. 
The position of the X-ray source (R.A., Dec. (2000) = 
13 h 29 m 54?17, +47°ll'30'/2) is fully consistent with the ra- 
dio position of the SN (~ 0"5 offset; Rupen et al. 1994). 
Astrometry of the M51 nucleus shows a similarly small off- 
set between the (2-8 keV band) X-ray and radio positions 
(~ 0'.'3). Exposure corrected source counts were extracted 
within a radius of 3 image pixels (90% encircled energy 
radius at 1 keV) and corrected for the background taken 
in an annulus with inner and outer radii of 3.5 and 9.5 
pixels, respectively (1 pixel corresponds to 0'.'49). 

A spectral model has to be assumed to convert source 
counts into energy fluxes. We adopted an effective (0.3- 

2 keV band) cooling function of A = 3 x 10 -23 ergs cm 3 s _1 
for an optically thin Raymond-Smith thermal plasma with 
a temperature of 10 7 K (Raymond, Cox & Smith 1976). 
Although the temperature is unknown, this value is con- 
sistent with theoretical expectations (Chevalier & Frans- 
son 1994), previous X-ray observations of other SNe and 
the detection of SN 19941 in the soft (0.3-2 keV) X-ray 
band where the peak of a 10 7 K spectrum is located. The 
equivalent count rate to (unabsorbed) flux conversion fac- 
tors are then 4 x 10 -11 (ergs cm~ 2 s _1 )/(counts s _1 ) and 

3 x 10~ 12 (ergs cm- 2 s- 1 )/(counts s" 1 ) for the ROSAT 
HRI and Chandra ACIS-S3, respectively, and a Galac- 
tic foreground column density of Ah = 1.3 x 10 20 cm~ 2 
(Dickey & Lockman 1990). The uncertainty of the conver- 
sion factors for optically thin thermal spectra with tem- 
peratures in the range 10 7 -TO 9 K is ~ 15%. Assuming a 
0.86 keV thermal bremsstrahlung spectrum (correspond- 
ing to a temperature of 10 7 K) instead of a Raymond- 
Smith thermal plasma increases the conversion factor by 
~ 7%. The Chandra results for the two individual obser- 
vations of SN 19941, together with previous ROSAT HRI 
results (Immler, Pietsch & Aschenbach 1998), are summa- 
rized in Table 1. 

3. DISCUSSION 

We used the Chandra and ROSAT data to construct a 
combined X-ray lightcurve of SN 19941, which is presented 
in Fig. 2. A best-fit single power-law rate of decline of 
/ x oc t~ s with index s — 0.9lg'2 (solid line, Fig. 2) is in- 
ferred for the measurements on days 82, 2,271 and 2,639. 
Assuming an initial exponential rise of the X-ray lumi- 
nosity after the outburst (at time to) of SN 19941 and a 
subsequent power-law decline with index s, we also param- 
eterized the X-ray evolution as / x oc (t — to)~ s x e~ r with 
t oc (t — to)~ & . This model has been successfully used to 
describe the time dependence of the radio emission of SNe 
(Weiler et al. 1996). The external absorption of the emis- 
sion is represented by the e~ T term ('optical depth') and 
the time-dependence of the optical depth is parameterized 
by the exponent [3. In the X-ray regime, the rise could 
represent either decreasing absorption by material along 
the line of sight to the hot gas or simply non-production 
of X-rays in the ROSAT band at early times. If we adopt 
this strictly heuristic description, we find s ~ 1 and (3^2 
(dashed line, Fig. 2). In case the ROSAT detection on day 
82 corresponds to a pre-maximum measurement, a steeper 



index of s ~ 1.5 is inferred (dotted line, Fig. 2). However, 
given the large error range of the flux measurements, s and 
(3 are not very well constrained. 

While a t~ x rate of decline has been observed for the 
Type lib SN 1999cm with Chandra (Pooley et al. 2002), 
faster rates of decline have been reported for the Ic SN 
1998bw (s = 1.4), inferred from the tentative hard (2- 
10 keV) band BeppoSax data (Pian et al. 1999), the Ic 
radio SN 1990B (s = 1.3) and the lb radio SNe 1983N 
(s = 1.6) and 1984L (s = 1.5; Van Dyk et al. 1993). 
By contrast, the long-term X-ray lightcurve of the lib 
SN 1993J is best described by a slow rate of decline with 
s = 0.27 (Immler, Aschenbach & Wang 2001). 

The X-ray emitting material could be either shocked 
SN ejecta (reverse shock) or shocked CSM (forward shock; 
Chevalier & Fransson 1994). The lack of time-dependent 
X-ray spectroscopy precludes a distinction but, as we shall 
see, the X-ray lightcurve is consistent with the latter model 
with a constant pre-SN mass loss rate (p csm oc r~ 2 ). Al- 
ternatively, the CSM might result from a non-conservative 
mass-transfer to the companion star, leading to the forma- 
tion of a flattened or disk-like H and He- rich shell (Nomoto 
et al. 1994). In this case a CSM profile flatter than 
Pcsm oc r -2 is expected due to the non-spherically sym- 
metric geometry of the mass-transfer. 

In the stellar wind scenario the continuum equation 
requires a mass- loss rate of M = Airr 2 p w (r) x v w (r) 
through a sphere of radius r. After the SN shock plows 
through the CSM, its density is p csm = 4p w (Fransson, 
Lundqvist & Chevalier 1996). The X-ray luminosity of 
the shock-heated CSM is L x = A(T)dVn 2 , where dV 
is the volume, n — p C sra/fn is the number density of 
the shocked CSM and m is the mean mass per particle 
(2.1 x 10~ 27 kg for a H+Hc plasma). Wc thus obtain 
L x = 4/(7rm 2 )A(T) x (M/v w ) 2 x (v s i) _1 . The observed 
X-ray rate of decline (t^ 1 ) is consistent with this descrip- 
tion if A(T), M/v w and v s are constant. We can hence 
use the observed X-ray luminosity at time t after the out- 
burst to measure the ratio M/v w assuming a constant shell 
expansion velocity v s . Given our observed X-ray luminosi- 
ties (Table 1) and assuming a shell expansion velocity of 
v s = 16, 500 km s _1 (Filippcnko et al. 1995) we derive a 
mass-loss rate of M ~ 1 x 1O~ 5 M yr _1 (u w /10 km s _1 ) 
consistent with all three detections and the ROSAT upper 
limit on day 1,368. 

It should be noted, however, that a single power-law 
model can only be reconciled with the data obtained on 
days 82-2,639, but not with the early ROSAT upper limit 
on day 52 (see Fig. 2). This is indicative that the above 
model might be incomplete for early epochs. The addi- 
tion of an exponential term describing the rise during the 
early phase of the emission leads to a fit that is in agree- 
ment with all measurements in case the ROSAT detection 
on day 82 was a post-maximum measurement (s ~ 1). By 
contrast, the model with a steeper rate of decline (s ~ 1.5), 
which postulates that the early ROSAT detection corre- 
sponds to a pre-maximum measurement, is in conflict with 
a Pcsm oc r~ 2 profile. 

Table 1 summarizes the CSM number density n for the 
different radii r = v s x t corresponding to the dates of the 
observations. As was already claimed based on the early 
ROSAT data (Immler, Pietsch & Aschenbach 1998), the 
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3cr upper limit for the CSM density at r = 2.0 x 10 17 cm 
(day 1,368) is lower than that at 1.2 x 10 16 cm (day 82). 
Our Chandra results at 3.2 and 3.8 x 10 17 cm (days 2,271 
and 2,639, respectively) clearly support these early indi- 
cations for a decreasing CSM density profile (see Table 1). 
The combined ROSAT and Chandra data of SN 19941 give 
a best-fit profile of p csm oc r _1,9±0 ' 1 . 

The only SN for which a CSM density profile has been 
constructed from X-ray measurements is SN 1993J (Type 
lib with a progenitor of ~ 15M Q ), based on long-term 
monitoring with ROSAT between six days and five years 
after the outburst (Immler, Aschenbach & Wang 2001). A 
comparison between the CSM profiles of the two SNe is 
presented in Fig. 3. Although the CSM density profile of 
SN 1993J is significantly flatter (p csnl oc r -16 ) than that 
of SN 19941, the two SNe have a similar CSM number den- 
sity of ~ 10 6 - 5 cm~ 3 at a radius of r ~ 5 x 10 15 cm from 



the site of the explosion (cf. Fig. 3). 

Given that both the X-ray rate of decline (L x oc t~ r ) 
and the CSM profile (p csm oc r- 19±01 ) of SN 19941 are 
not in conflict with what is expected for a constant stellar 
wind speed and a constant mass-loss rate of the progenitor 
(L x oc i -1 , p csm oc r~ 2 ), it is likely that the stellar wind of 
the massive progenitor dominates the CSM. 
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Table 1 



X-ray Properties of SN 19941 



Day < a ) 


Instrument 


Count Rate 






M W 


log n'°) 






(10~ 4 cts s- 1 ) 


(io- 15 ) 


(10 37 ) 


(10-S) 


(cm" 3 ) 


52 


ROSAT HRI 


< 3.1 


< 12.9 


< 10.9 


< 0.5 


< 5.9 


82 


ROSAT HRI 


5.2 ± 1.7 


21.6 ± 5.0 


18.2 ±4.2 


0.7 ±0.4 


5.7 


1,368 


ROSAT HRI 


< 8.4 


< 34.8 


< 29.4 


< 3.8 


< 4.0 


2,271 


Chandra ACIS-S 


5.5 ±2.2 


1.7 ±0.7 


1.4 ±0.5 


1.1 ±0.7 


3.0 


2,639 


Chandra ACIS-S 


4.6 ± 1.6 


1.4 ±0.5 


1.2 ±0.4 


1.1 ±0.7 


2.9 



(°) day after the outburst of SN 19941 (March 31, 1994; Chandler, Phillips & Rupen 1994) 

O 0.3-2 keV band flux in units of 10~ 15 ergs cm~ 2 s~ 1 for a 10 7 K thermal plasma spectrum and 
a Galactic foreground column density of TVjj = 1-3 X 10 20 cm~ 2 (Dickey & Lockman 1990) 

( c ' 0.3-2 kcV band luminosity in units of 10 37 ergs s _1 for an assumed distance of d = 8.4 Mpc 
(Feldmeier, Ciardullo & Jacoby 1997) 

' d ) pre-SN mass-loss rate in units of 10 — 5 Mq yr _1 

( e ) post-shock CSM number density 
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Fig. 1. — Chandra ACIS-S soft (0.3—2 keV) band X-ray image of the central region of M51. The image was adaptively smoothed to achieve 
a signal-to-noise ratio in the range 3-5 and is plotted in logarithmic greyscale. The position of SN 19941 is marked by a box and the nucleus 
of M51 is at the center of the image. 




Fig. 2. — Soft (0.3—2 keV) band X-ray lightcurve of SN 19941. The Chandra detections arc marked by filled boxes, open boxes indicate the 
early ROSAT HRI detection and 3cr upper limits. Error bars are ±lcr statistical errors. Time is given in days after the outburst. The solid 
line represents a power-law fit to the detections at days 82, 2,271 and 2,639. The dashed and dotted lines represent models including an e~ T 
term (see Section 3). 
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Fig. 3. — Circumstcllar matter density profile as a function of SN shell expansion radius. The solid line gives the best-fit CSM density 
profile of pesm tx r -1 ' 9 for the ROSAT (open boxes) and Chandra (filled boxes) measurements of SN 19941. The CSM density profile with 
Pcsm tx r -16 based on ROSAT observations of SN 1993J is drawn for comparison (dashed line and open triangles; Immler, Aschenbach & 
Wang 2001). 



